Draft version February 2, 2008 

Preprint typeset using I4TgX style emulateapj v. 10/10/03 



WHAT IS L*?: ANATOMY OF THE GALAXY LUMINOSITY FUNCTION 

ASANTHA COORAY 1 AND MlLOS MlLOSAVLJEVIC 2 
'Department of Physics and Astronomy, University of California, Irvine, CA 92617 
2 Theoretical Astrophysics, California Institute of Technology, Pasadena, CA 91125 
Draft version February 2, 2008 

ABSTRACT 

Using the empirical relations between the central galaxy luminosity and the halo mass, and between the total 
galaxy luminosity in a halo and the halo mass, we construct the galaxy luminosity function (LF). To the lu- 
minosity of the central galaxy in a halo of a given mass we assign log-normal scatter with a mean calibrated 
against the observations. In halos where the total galaxy luminosity exceeds that of the central galaxy, satellite 
galaxies are distributed as a power-law in luminosity. Combined with the halo mass function, this description 
reproduces the observed characteristics of the galaxy LF, including a shape consistent with the Schechter func- 
tion. When all galaxies are included, regardless of the environment or the Hubble type, the Schechter is 
the luminosity scale above which the central galaxy luminosity-halo mass relation flattens; L* corresponds to 
~ 10 13 M Q on the halo mass scale. In surveys where central galaxies in massive clusters are neglected, either 
by design or because of the cosmic variance, L + is simply the mean luminosity of central galaxies in halos at 
the upper end of the selected mass range. The smooth, exponential decay of the Schechter function toward high 
luminosities reflects the intrinsic scatter in the central galaxy luminosity-halo mass relation. In addition to the 
LF, the model successfully reproduces the observed dependence of galaxy clustering bias on luminosity. 
Subject headings: cosmology: observations — cosmology: theory — galaxies: clusters: general — galaxies: 
formation — galaxies: fundamental parameters 



1. INTRODUCTION 

The mass function of dark matter halos is routinely mea- 
sured in numerical simulations. If one assumes that the galaxy 
luminosity is proportional the halo mass, the abundance of 
galaxies at the faint and the bright ends o f the luminosity 
range is significantly below the expec ted ( Vale & Ostriker 
2004; van den Bosch, Yang, & Mo 2004). Semianalytic mod- 
els of galaxy formation, including those attempting to account 
for feedback and heating processes, generally do not explain 
the shape of the g alaxy luminosity function (LF; see, e.g., 
iBensonet all2003h. ' 

In lCoorav & M ilosavljevic1 J20051) . we studied the relation 
between the luminosity of the central galaxy and the mass 
of dark matter halo occupied by the galaxy (henceforth, the 
L c -M relation). We suggested that the flattening of this rela- 
tion in halos above M cr ; t ~ 1O 13 M0 is a consequence of the 
progressive decline in efficiency with which central galax- 
ies accrete satellites in the course of hierarchical merging. 
We demonstrated that to explain the L c -M relation, ongo- 
ing galaxy growth in massive halos via gas cooling and star 
formation need not be invoked. Therefore, we depart from 
the classical pic ture in which giant galaxies accrue mass con - 
tinuously (e.g., iRees & Ostrikerll 19771 iWhite & Reesll 19781) . 
and support the picture in which a bulk of the stellar mass is 
produced in smaller halos in which virial shoc ks do not raise 
the gas temperature to th e virial temperature (B innevlf2004t 
iDekel & Birnboim 2004). 

Here, we demonstrate that the luminosity function of galax- 
ies can be derived from two simple premises: (1) The lumi- 
nosities of central galaxies in halos of a given mass possess 
a log-normal intrinsic scatter with a mean and a dispersion 
identified with measured values, and (2) When the total lumi- 
nosity of galaxies inside a halo exceeds that of the central 
galaxy, the luminosities of satellite galaxies are distributed 
as a power-law. These two assumptions yield a LF of the 



lSchechtert ( fT97^) type, $(L) oc L a e ' *, where a is the slope 
of the LF at the faint end, and is usually thought of as a 
characteristic luminosity scale. 

Our construction of the LF is compatible in spirit with the 
ab initio synthesis of the conditional stellar mass function in 
Zhen g"etall J2004I) . and differs fro m the modeling o f condi- 
tional luminosity function (CLF) in lYang et alj j2005), where 
the Schechter form was assumed a priori. Log-normal scat- 
ter of galaxy luminosities in low-m ass halos was invoked by 
Yang, Mo, & van den Bosch (2003) to explain the scatter in 
the Tully-Fisher (TF) relation JTullv & Piercdl2000i) . Here, 
we show that the log-normal scatter is crucial for explain- 
ing the turnover of the Schech ter LF at the brieht-end . Our 
approach is opposite to that of iVale & OstrikerT (f2004). who 
used the LF to extract information about the L c -M relation. 
We instead use the observationally determined L c -M relation 
to reconstruct the LF. We also calculate the CLF, and illustrate 
how the LF is built from the CLF. 

In § |2j we describe the construction of the galaxy LF. In 
§|5] we compare our LF with the observed LFs of cluster and 
field galaxies. We adopt the curre nt concordance cosm ologi- 
cal model consistent with WMAP ( Spergel et al. 2003) with a 
Hubble constant of h = 0.7 unless stated otherwise. 

2. MODEL LUMINOSITY FUNCTION 

The CLF, denoted by $(L|M), is the average number of 
galaxies with luminosities between L and L + dL that reside 
in halos of mass M. We separate the CLF into terms associ- 
ated with central and satellite galaxies, $(L|M) = $ C (L|M) + 
<P S (L|M), where the central galaxy CLF is a log-normal dis- 
tribution in luminosity with a mean L C (M) and dispersion 
ln(10)S 

* e( L|M)= * (M) ^L^oiL/umn (1) 

1 V2^1n(10)EL P \ 2£ 2 J 

and the satellite CLF is a power law $ S (L|M) = A(M)L 7 , 
where $(M) and A(M) are normalization factors. The sepa- 
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FIG. 1. — (a) C entral galaxy luminosity as a function of the halo mass. The data are model fits to the SDSS galaxy-mass correlation function 
i Cooray & Milosavlievic 2005; red points), the same masses estimated by Yang et al. (2003; blue squares), a direct measurement of galaxy luminosity and 
halo mass for a sample of galaxy groups and clusters from Lin & Mohr (2004; cyan points). To extend the relation below the luminosities considered by 
Coorav & Milosavlievic 1 2005 ), we use the low luminosity data points from Vale & Ostriker (2004; magenta squares). We also show the total luminosity of 
galaxy groups and clusters based on data of Lin et al. (2004; small squares), (b) The LF of galaxies (thick solid line). We show the contributions to CLF from 
central galaxies dashed red lines. From left to right, the CLFs are shown for dark matter halos with masses separated into six decade intervals spanning the 
range (10 10 — 10 16 )Mq with E = 0.25. The total LF is shown as a solid black line. The blue line is the total LF related to satellite galaxies in halo masses 
above ~ 10 1 *Mq with 7 = — 1. Central galaxies dominate the total LF at any luminosity. For comparison, we also show the LF obtained assuming a constant 
mass-to-light ratio M/L = 100 (dot-dashed line). The total LF can be approximated by a Schechter function with a fa —1.3 and L* fa 3 X 1O"L0, where a is 
measured at L = 1O 8 Z,0 . The value of L k shifts to lower luminosities when one ignores central galaxies in groups and clusters. 



ration into central and satellite galaxies is motivated by the 
halo model for galaxy statistics (Cooray & Sheth 2002|). 

The normalization $(M) of the central galaxy CLF is fixed 
such that J <&{L\M)LdL equals the average total luminosity 
of galaxies L tot {M) in a halo of mass M. The normalization 
A(M) of the satellite CLF can be obtained by defining L S (M) = 
Ltot(M) -L C (M) and requiring that L S (M) = ®s(L\M)LdL, 
where the minimum luminosity of a satellite is L m j n , while we 
generally set L max = L c . We find below from an analysis the 
cluster LF that a more appropriate value for the maximal lu- 
minosity of satellites is between L c /2 and L c /3. Our recon- 
struction is independent of the exact value assumed for L m j n , 
as long as it lies in the range (10 6 - 1O 8 )L . Note that L tot (M) 
equals L C (M) for M < 10 n M Q , and thus L S (M) = on these 
scales. The remaining two free parameters are 7 and S; we 
discuss the determination of these parameters in §0 

We adopt the Sheth & Tormen (1999; ST) mass function 
dn/cLM for dark matter h alos as it approximates n umerical 
simulations better than the Press & Schechterl (119741) function 
(see, e.g.. l.Tenkins et al.l2001l) . The usual, unconditional LF is 
then given by 

HD = [ HL\M)^dM. (2) 

The CLF represents galaxy statistics better than the LF when 
wide-field data sets are available in which redshiftsare mea- 
sured for tens of thousands of galaxies (Yang et al. 2005). 

For L C (M), we make use of the observed central galaxy lu- 
minosity measurements over six orders of magnitude in l umi- 
nosity (see Fig. [^). Following Vale & Ostriker (2004), we 
employ a fitting function of the form 

(MlMx) a 

L(M) = Lq v ' . (3) 

[b + (M Mi) cd ] l / d V ; 



For central luminosities, the parameters are Lq = 4.4 x IQ 11 L Q , 
Mi = 10 n M Q , a = 4.0, b = 0.9, c = 3.85, and d = 0.1 . 
These values are different than in iVale & Ostrikerl ( 120041) : 
their fit was constructed from the fey-band luminosities of 
2dF, while we use the K-bwA luminosities. For L > 
5 x 10 10 Lf7), the relation is co mpatible with that derived in 
iCoorav & M ilosavlievic (2005). For total luminosities, we 
also use the fitting formula in equation 0, but with c = 3.65. 
At the massive end, the total luminosity can alternatively 
be described as a power-law. The constructed LF does not 
change if power-law behavior is enforced there. This is be- 
cause the total LF (including field, group, and cluster galax- 
ies) is dominated by central galaxies on any scale. The overall 
shape of the LF is thus sensitive to the shape of the L c -M re- 
lation. 

3. RESULTS AND DISCUSSION 

In Figure ^i, we show the constructed LF in which the 
galaxies in halos in the mass range (1O 1O -1O 16 )M have been 
taken into account. We also plot the contributions to CLF 
from central galaxies and the contribution to LF from satellite 
galaxies. The central galaxies dominate t he LF on any lumi- 
nosity scale, similar to the conclusion bv lZheng et alJ ( 120041) 
that central galaxies dominate the stellar mass function on any 
mass scale. The fractional contribution of satellites to LF de- 
pends on the luminosity and has a maximum slightly below 
~ 3 x 10 u Lq. We now explore the origin of the slope a 
of the LF at the faint end, and then discuss the meaning of the 
characteristic scale above which LF decreases steeply with 
increasing luminosity. 

The faint-end slope of the LF is related to the low-mass 
slopes of the halo mass function and the L c -M relation. At 
low halo masses, the mass function dn/dM oc M~ 2 a(M)~^ , 
where (3 = 1 in the Press-Schechter mass function and j3 w 0.4 
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FIG. 2. — (a) Galaxy cluster LF for 7 = — 1, where luminous central galaxies are ignored in the construction (blue lines), for varying maximum luminosity of 
a satellite L raax = (1, 1/2, 1/ 3,1/5)L C (blue lines). The cu rves are arbitrarily normalized at the faint end of luminosity as the absolute normalization of 4>(k) is 
unknown. The data are from Lin, Mohr. & Stanford 1 2004). For reference, we also show the contribution to the LF by central galaxies (red dashed line). The total 
LF clearly departs from the Schechter form, (b) The A'-band field galaxy LF (black lines); we show the galaxy LF with maximum halo mass of (0. 1 , 1) X 10 13 Mq 
and scatter S = 0.23. The dot-dashed line is the LF with all halos included. The thin red lines with show the variation of the LF when the width of the log-normal 
distribution, E, varies to 0.15 and 0.3, assuming a maximum mass of 1O 13 M0. Note that S PS 0.23 describes the data best. The data are from Huan g et alj 120031) . 
We also shows the Schechter function fit to the data in Huang et al. 1 2003) with Mk* = -23.7, a = -1.38, and = 0.013 Mpc~ 3 (dashed blue line). Note that the 
Hubble constant h = 1 in panel (b). 
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FIG. 3. — Galaxy bias as a function of luminosity calculated from con- 
ditional luminosity functions (black line); we include the separate contribu- 
tions from central galaxi es (red line ) and satellites (blue line). Also shown 
are the SDSS of Zehavi et al. 1 2004), and a fit to the galaxy bias in 2dF data 
of Norberg et al. 1 2002) (green dashed line). 



in the ST mass function. Here, cr(M) oc M" ( " +3) / 6 is the den- 
sity variance on scales M, while n is the slope of the matter 
power spectrum, P(k) oc k". Since central galaxies dominate 
the LF, we have $(L) ~ / $ c (L\M)(dn/dM)dM. If L oc M", 
we can write $(L) - J L'~ 1 ~^ v+l3{ " +3) / 6ri S(L' - L)dL' , where 
we have ignored the scatter in the L c -M relation by setting 
£ — > and thus replacing the log-normal distribution repre- 
senting the dispersion of central galaxy luminosities with a 
(5-function at L c ; the faint-end slope of the LF is insensitive to 
E. The faint-end LF then scales as $(L) oc L-i-i/V«"+3)M 



An examination of Figure^ shows that the faint-end slope 
of the LF is defined by galaxies in dark matter halos with 
masses below ~ 10 n M Q . At these mass scales, 2 < 77 < 4, 
where the upper limit corresponds to asymptotically low mass 
halos. Since -2 < n < -1, we have f3{n + 3)/6r] <C 1, and 
thus $(L) ~ L a with a > -1.5. We expect a < -1.25, un- 
less the L c -M relation is steeper than 77 = 4 at the low-mass 
end. We expect that the slope of at the faint end is 
independent of the color selection and bandpass, as these 
choices only affect the normalization of the mass-to-light 
scaling. In the case of the /T-band gala xy LF, the measu red 
faint-end slopes r ange from a = - 1 .39 (Huang et al. 2003) to 
a = -0.93 ± 0.04 (ICole et al.l2001l) . At other wavelengths, the 
best-me asured galaxy LF ha s a = —1.21 ±0.03 in the bj band 
of 2dF (Norber g et alJl2002l) and the flatter a = -1 .05 ± 0.01 
in the r'-band of SDSS JBlanton et al.l2003T) . When the same 
luminosity func tion was e xtended to fainter magnitudes, how- 
ever, a = —1.3 (Blanton et al. 2004), consistent with our ex- 
pectations. 

Provided that the L c -M relation is well known on all scales, 
the LF of galaxies subject to a particular selection criterion 
based on galaxy type or environment (e.g., field galaxies, clus- 
ter galaxies, etc.) can be estimated directly by inserting the 
probability of selection in the integral in equation (|2j- Vari- 
ation between selection criteria may be responsible for the 
observed variation in a. Selection criteria also seem to be re- 
sponsible for the measurements with a > — 1, which naively 
seem discrepant with the above predictions. If, e.g., a survey 
is biased toward ignoring the galaxies in low mass halos, one 
may indeed measure a > — 1. 

It is unclear whether the faint-e nd slope of the LF changes 
with redshift. Dror v et al] (12003). who studied the redshift 
evolution of the LF in the /^-band, assumed a priori that a is 
independent of redshift. Since the slope of the mass function 
at low masses is independent of redshift, any variation of a 
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with redshift would imply that the L c -M relation must also 
vary with redshift. The expected correlation between a and 
the slope of the L Q -M relation can be tested directly if a galaxy 
survey, e.g., atz~ 1, can be combined with weaklensing mass 
estimates for galaxies in the same sample. 

When an unbiased galaxy sample is used, L± can be iden- 
tified with the scale above which the luminosity scatter in 
the L c -M relation dominates over the increase of luminos- 
ity with mass, i.e., dlnLjdlnM ~ ln(10)£. With £ - 0.25, 
this yields M t ~2x 10 13 M Q and U = L C (M„) ~3x 10 n L Q . 
Above this value, the further increase in mean galaxy lumi- 
nosity with halo mass is smaller than the scatter, and the cor- 
responding CLFs overlap. Had the L c -M relation continued to 
increase at the rate dlnL c /dlnM > ln(10)£, the CLFs would 
have been disjoint and the Schechter function would not have 
exhibited the usual exponential cut-off. 

Satellite galaxies do not contribute significantly to the over- 
all LF. Nevertheless, galaxy cluster studies routinely mea- 
sure the LF of satellite galaxies in clusters. In Figure 
Ek, we plot the LF of sate llite galaxies in clusters from 
iLin. Mohr. & Stanford ll2004l) . Since the measurement is in 
the Zf-band and is based on the same data is in Figure^, a 
direct comparison with our model is possible. In Figure^, 
we plot the LF with and without central galaxies, and also 
the LF of central galaxies alone. In this calculation we as- 
sumed that the cluster mass lies in the range (1O 13 -1O 15 )M 
with a dispersion of £ = 0.25. Our construction re produces 
the observational fact llLin. Mohr. & Stanforcll2004l) that the 
cluster LF, inclusive the central gal axies, departs sigrrificantly 
from the Schechter form. Similarly, Trentham & Tullv (2002) 
decomposed the total LF of nearby clusters and groups into 
a log-normal and a Schechter function-like component. The 
ILin. Mohr. & Stanford! J2004) data suggest that the power-law 
slope of the satellite distribution is 7 w — 1 ±0.2, while the 
maximum luminosity L max of the satellite distribution is in the 
range (1/3-1 12)L C . 

In Figure |2p, we compare the field galaxy LF of 
iHuang et alJ (|2003) with the predictions of our model. The 
definition of field galaxies here includes all galaxies regard- 
less of the environment. Nevertheless, their selection of 
galaxies may still somehow have been biased or affected 
by the cosmic variance (the random fluctuation in the num- 
ber of groups and clus ters in a survey of limited volume). 
The Huang etal] (|2003) data are best described with a max- 
imum mass of ~ 1O 13 M . This corresponding is fac- 
tor of 2 below the value expected based on the construc- 



tion described in § [2] The discrepancy could perhaps be as- 
cribed to the cosmic variance in t he number of clusters in 
the sample of IHuang etail II20031 We vary the width £ 
of the Gaussian that describes the dispersion in luminosity 
at fixed halo mass. To match the exponential behavior of 
the LF, £ ~ 0.23 is appropriate, which is perfectly consis- 
tent with the factor of 2 vertical scatter in the L c -M relation. 
This scatter is somewhat higher tha n £ ~ 0.17 estimated by 
Yang. Mo. & van den Bosch (2003) for the bright-end of the 
TF relation (mass scales of 1 13 Mf^); their estimate of the 
scatter was based on the data of Tullv & Pierce ( 2000). 

As an example of a practical application of our re- 
constructed CLF, in Figure [5] we compare our pre- 
diction of the luminosity-dependen t galaxy bias with 
the m easurements f rom S PSS ( Ze havi et alJ 120041) and 
2dF (Norberg etail 120021) . The bias equals b(L) = 
<J>(L) _1 J b{M)${L\M)(dn / dM)dM, w here b(M) is the halo 
bias b ased on the ST mass function ( Shet h. Mo. & TormenI 
12001 . The bias increases monotonically from below unity 
to above unity as L/L* ranges from 0.01 to 10. This behavior 
is reproduced in our model which shows that the average bias 
is also dominated by central galaxies at a given luminosity. 

Semi-analytic studies h ave, thus far, not suc ceeded in repro- 
ducing the observed LF (Benson et al. 2003). Here, we have 
presented a simple model that exposes the basic elements that 
shape the LF. The faint-end slope of the LF has received at- 
tention recentl y because it differs from the slope of the halo 
mass function (Ben son et alJ l2002). If the scaling of the av- 
erage galaxy luminosity with the halo mass at luminosities 
below can be explained, LF can be directly recovered as 
we demonstrate in § [2] The shape of the LF at L > is 
then governed by the dissipationless mergi ng model for lu - 
minosity growth presented in Coorav & Milosavlievic ( 2005 ). 
There, we argued that feedback processes must prevent con- 
tinued star formation on all luminosity scales > L*. We also 
note that our construction may also make it easier to reliably 
measure cos mologic al parameters from the galaxy LF, as first 
attempted in Seliak (2002). 
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